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Abstract 

Recent experimental results on exclusive semi-tauonic B meson decays, B — > D^*'rv, showing 
sizable deviations from the standard model prediction, suggest a new physics in which the structure 
of the relevant weak charged interaction may differ from that of the standard model. We study 
the exclusive semi-tauonic B decays in a model-independent manner using the most general set 
of four-Fermi interactions in order to clarify possible structures of the charged current in new 
physics. It turns out that correlations among observables including tau and D* polarizations and 
q 2 distributions are useful to distinguish possible new physics operators. Further, we investigate 
some interesting models to exhibit the advantage of our model-independent analysis. As a result, 
we find that two Higgs doublet models without tree-level flavor changing neutral currents (FCNC) 
and the minimal supersymmetric standard model with i?-parity violation are unlikely to explain 
the present experimental data, while two Higgs doublet models with FCNC and a leptoquark model 
are consistent with the data. 
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I. INTRODUCTION 



The V — A structure of the weak interaction was established in the study of nuclear (3 
decays for the quarks and lepton in the first generation |T]. Since then, several charged 
current interactions in the second and third generations have been investigated and turned 
out to be described by similar V — A interactions with a reasonable level of accuracy. 

Among the charged current processes observed so far, pure- and semi- tauonic B meson 
decays, B — )■ rv [21 [3] and B — y D^*'tv [IHH] are thought-provoking in the sense that they 
contain both the quark (b) and the lepton (r) in the third generation. Their large masses 
suggest a close connection with the electroweak symmetry breaking (EWSB), which is still 
not fully understood even after the discovery of a di-gamma resonance [9j [TO] ■ 

The EWSB sector in a candidate of new physics beyond the standard model (SM) often 
has a different structure from that of the SM. A typical example is the two Higgs doublet 
model (2HDM) of type II [11], which has the Higgs sector identical to that of the minimal 
supersymmetric standard model (MSSM) [12] at the tree level. A pair of charged Higgs 
bosons appears in this class of models, and its couplings to fermions are proportional to the 
involved fermion masses and further enhanced if the ratio of vacuum expectation values, 
tan/3, is large. Several authors have studied effects of these enhanced couplings, which 
modify the V — A structure, in the pure- and semi- tauonic B decays [T3H26] . 

From the experimental point of view, these B decay processes are rather difficult to 
identify because of two or more missing neutrinos in the final states. With the help of 
large statistics and low backgrounds, however, they are good targets of e + e~ B factory 
experiments. The observed branching fraction of B — y tv is O(10~ 4 ) and that of B — y D^*'tp 
is O(10~ 2 ). Comparing the pure- and semi- tauonic B decays, the latter provides a wide 
variety of observables besides the larger branching fraction, such as decay distributions 
[HI [20H2H I2ZHS] and polarizations [151 1251 123 12H] • Hence, the semi-tauonic processes allow 
us to investigate the structure of the relevant charged current interaction, and we concentrate 
on B — y D^tv in this work. 

Another advantage of B — y D^*'rv is that theoretical and experimental uncertainties 
tend to be reduced by taking the ratio of branching fraction to the semileptonic decays, 
B — y D(*>£&, where £ denotes e or /i. These decay modes including lighter charged leptons 
are supposed to be less sensitive to the EWSB sector and assumed to be described by the 
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SM in the present work. 

The ratios of the branching fractions are defined by 

R{D = B(B -)• DW-P t ) ' (1) 
and their experimental values are summarized as 

R(D) = 0.42 ± 0.06, R(D*) = 0.34 ± 0.03, (2) 

where we combine the results of BaBar [5] and Belle [6HH] assuming the gaussian distribution. 
These measurements seem to be larger than the SM predictions, 

R(D) = 0.305 ± 0.012, R(D*) = 0.252 ± 0.004, (3) 

and the SM is disfavored at 3.5cr. 

Interestingly, as is shown in Ref. [5], one of R(D) and R(D*) is made consistent with the 
corresponding experimental value in the 2HDM of type II by adjusting tan/3/m#±, but it 
is unlikely to make both of them agree with the experimental results. Thus, the 2HDM of 
type II is also disfavored. 

In the present work, we first provide a model-independent framework for the semi-tauonic 
B decays in order to clarify the type of new physics that is required to explain the current 
experimental results. Then, we apply the framework to analyzing several models. The 
rest of this paper is organized as follows: In Sec. [TTJ we present the most general effective 
Lagrangian of b — > crv and the resultant helicity amplitudes of B — > D^*>tv. Experimental 



constraints and theoretical predictions are given in Sec. Ill We discuss how to distinguish 



new physics contributions that appear in the effective Lagrangian in Sec. IV We analyze 
2HDMs, MSSM with i?-parity violation and a leptoquark model in Sec. |VJ Section VI is 
devoted to our conclusions. 



II. EFFECTIVE LAGRANGIAN AND HELICITY AMPLITUDES 

A. Effective Lagrangian for b — > crv 

The semi-tauonic bottom quark decay, b — > crv, is described by the four-Fermi interaction 
of charged left-handed currents in the SM. Other four-Fermi operators might be induced in 
the presence of new physics. 
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The effective Lagrangian that contains all conceivable four-Fermi operators is written as 
- £ cff = 2V2G F V cb [fa + C l Vl )O l Vi + C l V2 O l V2 + C l Si O l Si + C l S2 O l S2 + C l T O l T ] , (4) 

l=e,/i,T 

where the four-Fermi operators are defined by 

O l Vl = CL'fbLTLl^VU, (5) 

O l V2 = c R ^b R f Llfl u Ll , (6) 

1 51 = c L b R f R u u , (7) 

1 52 = c R b L f R u Ll} (8) 
O l T = c R a^b L f R a^u Ll , (9) 

and C l x (X = Vi j2 , Si : 2,T) denotes the Wilson coefficient of O l x . Here we assume that the 
light neutrinos are left-handed. 1 The neutrino flavor is specified by /, and we take all cases 
of / = e, fi and r into account in the contributions of new physics. Since the neutrino flavor 
is not observed in the experiments of bottom decays, the neutrino mixing does not affect the 
following argument provided that the Pontecorvo-Maki-Nakagawa-Sakata matrix is unitary. 
The SM contribution is expressed by the term of Si T in Eq. Q. We note that the tensor 
operator with the opposite set of quark chiralities identically vanishes: cia^ v b R f R a^ v uu = 0. 



B. Helicity amplitudes 

The helicity amplitude of B — > Mtv (M = D,D*) is written as 

Mf T ' XM = 5 lT M^ Am + Mylf M + Mylf M + M x s [f M + Mg T 2 f M + M x T T f M , (10) 

where A T is the helicity of the tau lepton, X M = s indicates the amplitude of B — > Dtu, and 
X M = ±1,0 denotes the D* helicity defined in the B rest frame . The amplitude M X l^ M 
represents the SM contribution and other terms in the right-hand side stand for new physics 
contributions corresponding to the operators in Eqs. (Jsj) — (Joj) - The SM amplitude is given by 

m 



c 

m Xt,Xm 



l SM 



^V^VxH^L^, (11) 



1 Possibilities to introduce a light right-handed neutrino are discussed in Refs. 
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and the amplitudes that represent new physics contributions take the following forms: 



M X y % f M = -^Vc b C l Vi £ VxH^ x L* :i (i = 1,2), (12) 
M X s \i M = -^V ch C l Si H^ -If (i = 1,2), (13) 



where if's and L's are the hadronic and leptonic amplitudes respectively as defined below; 
A, A' = ±,0,s denote the virtual vector boson helicity. The metric factor r]\ is given by 
ri± = 1 and rj s = —1. We treat the contraction of the Lorentz indices in O l T as if two heavy 
vector bosons are exchanged. 

The leptonic amplitudes, L x T t , Lf T and L x T y z are defined by 

L^(q 2 lC os6 T ) = e M (A)(r(p T ,A r )^(^)|r7 At (l-75)^|0), (15) 
L^(q 2 ,cos9 T ) = (T(p r ,A r )^)l^(l-75H0>, (16) 
L^ ^q 2 , cos 9 T ) = -^ t (A)e,(A / )(r(p r ,A T )^(^)|r(T^(l- 75 )^|0), (17) 

where e M (A) represents the polarization vector of the virtual vector boson, = p B — p%{= 
p^ + Py) is the momentum transfer, and 9 T denotes the angle between the momentum of 
the tau lepton and that of the M meson in the rest frame of the leptonic system, to which 
we refer as the q rest frame [2S]- The r helicity A T is also defined in the q rest frame. The 
explicit formulae of the leptonic amplitudes are relegated to Appendix [Aj 
The hadronic amplitudes, Hy' x , H S M and H xx , are defined by 





= e;(X)(M(p M ,e(X M )) \crf{l - >f)b\B(p B )) , 


(18) 




= e;(\)(M(p M , e (A M )) |c 7 m (1 + l")b\B(p B )) , 


(19) 


H^(q 2 ) 


= (M(p M ,e(AM))|c(l + 7 5 )&|5( P B)), 


(20) 




= (M(p M ,e(X M ))\c(l - j 5 )b\B(p B )) , 


(21) 




= ie;(X)e:(X')(M(p M ,e(X M )) \ca^{l - ^)b\B{p B )), 


(22) 



where e(Xu) denotes the polarization vector of D* for Xm = ±1,0. The relations H v x = 
Hy 2 x and = H S S2 hold because of the parity conservation in the strong interaction. 
Similarly, we find H X S M = -H* M for X M = ±1,0. 
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The hadronic amplitudes of the vector type operators for B — > Dtv defined in Eqs. (18) 



and (19) are represented as 

H^ ± (q 2 ) = H^ ± (q 2 ) = 0, 
H s Vifl {q 2 ) = ^ 2i0 (g 2 )=m BV ^(l + r 

H^ s {q 2 ) = H^ s (q 2 ) = m B ^{\ - r 
and those for B — > D*tu are 



Vw 2 -1 

Vv 2 ( w ) 



w. 



Si(w) 



Hv u ±(<l 2 ) = -H^ T {q 2 )=m B ^A 1 {w) 
m B y/w 2 — 1 



q 2 (w 



w + 1 =F Vw 2 - lR^w) 



--Ax{w) [-w + r + (w - l)R 2 (w)] 



(23) 
(24) 

(25) 

(26) 
(27) 



q 2 (w) 



Ai(iu) [-2r(w + 1) + (1 - r 2 )R 2 (w) - q 2 (w)R 3 (w)] , (28) 



others = 0, 



(29) 



where q 2 {w) = 1 + r 2 — 2rw, r = mm/mB, and w = p B ■ PAi/(m B mM) is the velocity transfer 
in B — > Mtv. The form factors, V\(w), Si(w), Ai(w) and Ri(w) are defined in Appendix [b| 
The amplitudes of the scalar type operators are expressed in terms of vector form factors 
by applying the equations of motion of the quark fields: 

HIM") = HI 2 (q 2 ) = m B ^(w + l)S 1 (w), 
HW) = Hi 2 (q 2 ) = 0, 



(30) 
(31) 



him 2 : 



2VF(l + r) MW) t~ 2r{w + 1} + (1 " r2)i?2(w) " ^H^M] 
Similarly, for the tensor operator, we obtain 



,/ ?/ ,2 _ 1 

^ X ^ r (l + r) 2 V 1 (w) + 2r(w + l)S 1 (w)} 



H s + _{q 2 ) = Hl{q 2 )=m B ^ 



H± {q 2 ) = ±H± s (q 2 ) 
m B ^/r 



q 2 (w) 



- Ai{w) ±(1 - r)(w + 1) + (1 + r)Vw 2 - lR^w) 

VQ 2 ( W ) L 
HUf) = H° 0s (q 2 ) 

= rnsVl + r){w + i) + (w 2 _ i )jR3H ] 

1 + r L J 

H^,{q 2 ) = -H X ^(q 2 ), 



(32) 

(33) 
(34) 

(35) 
(36) 
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and others = 0. A detailed derivation of these hadronic amplitudes is found in Appendix |B} 



C. Form factors 

The form factor V\(w) is extracted from the experimental data on B — > D£u, where 
£ — e, p, provided that the decay process is described by the SM. Employing the following 
ansatz [36J, 

Vl(w) = Vi(l) [l - Splz + (51.pl - 10.)z 2 - (252. pj - 84.)^ 3 ] , (37) 

where z = (y/w~+l — \/2)/(y/w + 1 + V%), the heavy flavor averaging group (HFAG) deter- 
mines the slope parameter p\ as p\ = 1.186 ± 0.055 [37] . 

The form factor Si(w) does not contribute to the semileptonic B decay into a massless 
charged lepton. However it reduces to the Isgur-Wise function j3H] in the heavy quark limit 
as well as Vi(w). Accordingly we parametrize it as 

S 1 (w) = [l + A(w)]V 1 (w), (38) 

where A(w) denotes the QCD and 1/rriQ corrections. Following Refs. [THJ, ES] , we estimate 
A(w) and give an approximate formula [25] : 

A(iu) = -0.019 + 0.041(w - 1) - 0.015(w - l) 2 . (39) 

As for B — > D*£u, three form factors, A\(w) and Ri^{w), contribute. They are extracted 
from experimental data using the following parametrizations: 

A x (w) = A X {1) [I - 8p 2 Ai z + (53/4 - 15 )^ - ( 231 PX ~ 91)^ 3 ] , (40) 
Ri(w) = Ri{l) - 0.12(w - 1) + 0.05(w - l) 2 , (41) 
R 2 (w) = R 2 {1) + 0.11(w - 1) - 0.06(w - l) 2 . (42) 



Eq. (40) is given in Ref. [36] as well as Vi(w). The w dependence of Ri^ is estimated by the 
heavy quark effective theory [39], while -Ri(l) and ^(1) are left as fitting parameters. The 
HFAG determines these parameters as follows [37]: p\ = 1.207±0.026, R X {1) = 1.403±0.033 
and ^2(1) = 0.854 ± 0.020. The form factor Rs(w) appears only in B — > D*tv and is 
estimated by the heavy quark effective theory as [36J 

R 3 (w) = l+A 3 (w), A 3 (w) = 0.22 -0.052(w- 1) + 0.026(w - l) 2 . (43) 
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III. EXPERIMENTAL CONSTRAINTS AND THEORETICAL PREDICTIONS 



A. Observables and experimental constraints 

There are several measurable quantities affected by new physics operators in B — > D^*'tv. 
We consider the decay rate and the tau longitudinal polarization in B — > Dtv and define 
the following quantities: 

r+( C) + r-( C ) = YHD) - r-(D) 

V ; T(B — >■ Dlv) ' V ; V+{D)+T-(D)' V ; 



where r ± (D) represents the decay rate with A r = ±1/2 as defined in Eq. (C3) and T(B — > 
Dlv) is the total decay rate of B — > Dtv. The uncertainties that come from V c b and the 
normalization factor Vi(l) vanish in the above quantities. For B — > D*tv we also define 
R(D*) and P T (D*) in the same fashion. In addition we introduce the D* polarization as 

T{Dj) 

D * T(D* L ) + r(£>* ) ' 1 ' 

where -D^(t) represents the longitudinally (transversely) polarized D* and T(D* LT ) are de- 
fined in Eq. (C4). The uncertainties in R(D*), P T (D*) and Pu* due to V c b and Ai(l) also 
disappear. The new physics operators are expected to affect these observables in various 
ways. Thus it is important to study them at the same time in order to distinguish the 
underlying new physics. 

Besides the above integrated quantities, q 2 distributions are potentially sensitive to new 
physics. As we will illustrate below, the q 2 distribution of B — > Dtv decay rate is helpful 
in discriminating between two scalar operators. 

In the following model-independent analysis, we assume that one of the new physics 
operators in Eq. (|4]) is dominant except the SM contribution. This assumption allows us 
to determine the dominant Wilson coefficient from the experimental results of R{D) and 
R(D*), and to predict other observables. A situation beyond this assumption is discussed 
in Sec. [V] 

In Figs. [I] and [2J we show new physics effects on R(D^), P T (D^) and Pjy*. The horizon- 
tal axis is \Cx\ and three cases of Sx = 0, tt/2, and ir are shown for illustration, where Sx is 
the complex phase of C x . Effects of O e x ^ are the same as those of O x with 5x = tt/2 because 
the new physics contributions do not interfere with the SM amplitudes in these cases. The 
width of each prediction indicates uncertainties due to the form factors. We include ±100% 
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FIG. 1: Predictions on the branching ratios as functions of the absolute value of Wilson coefficient 
\C T X \ for X = Vi t 2-i S\^,T. The predictions of new physics effects for the operators O x ^ are given by 
the lines for Sx = vr/2 in these graphs. The light blue horizontal bands represent the experimental 
values. 

errors in the overall magnitudes of A(w) and As(w) as uncertainties in addition to the ranges 
of p\,p\ v -Ri(l) and i?2(l)- The horizontal bands with dashed boundaries in Fig. [ljrepresent 
the experimental values given in Eq. (|2]). From these results, we find that the sensitivity to 
the magnitude of the Wilson coefficient varies depending on each operator. We note that 
the theoretical uncertainties are sufficiently smaller than the present experimental accuracy. 
Therefore, we use the central values of the theoretical predictions in the rest of this work 
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FIG. 2: Predictions on P T (D),P T (D*) and P D * as functions of \C T X \ for X = Vi >2 ,S lj2 ,T. The 
predictions of new physics effects for the operators O e x are given by the lines for 5x = tt/2. 

unless otherwise stated. 

In Fig. [3j the allowed region of the complex Wilson coefficient C T X is shown for each 
operator O x . The vector operators O v describe the current experimental results [32]. The 
operator O t s is unlikely while O t S2 is favored as is already pointed out in Refs. [32], HQ]. In 
addition, we find that the tensor operator O r T reasonably explains the current data. We can 
read the allowed regions of \C X ^\ from that of the pure imaginary C x , since neither C e x nor 
the pure imaginary C r x interferes with the SM contribution as is mentioned above. Thus, 
we find no allowed region of the Wilson coefficient within 99% confidence level (C.L.) for 
Og^ nor Og^. The operators Oy, Oy^ and Ojf 1 are able to explain the present data. 



B. Correlation between R(D) and R(D*) 

Since each new physics operator contributes to R{D) and R(D*) in different ways as seen 
in Fig. [TJ it is useful to examine the correlation between R(D) and R(D*) for the sake of 
discrimination of new physics operators. In the left panel of Fig. |1| we show R(D) and 
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FIG. 3: Constraints on the Wilson coefficient C x . The allowed regions within 90% (blue), 
95%(cyan) and 99%(light blue) C.L. are shown. Allowed regions for C^H are obtained by those 



for pure imaginary C x . The leptoquark panel is mentioned in Sec. VC 



R(D*) in the presence of the new physics operators for O t x (X = Vi^Sip and T). The 
shaded regions are predicted by the indicated operators. As is known [15], R{D) is more 
sensitive to the scalar type operators O t s and O t S2 than R(D*). This is due to the angular 
momentum conservation that gives an extra suppression factor for the vector meson. On the 
other hand, we find that the tensor type operator OTp exhibits the opposite behavior, that 
is, R(D*) is more sensitive to Oj> than R(D). The vector type operator O v gives a unique 
relation between R(D) and R(D*), since it is just the SM operator and changes only the 
overall factor. The vector type operator Oy 2 that contains the right-handed quark current 
covers a wide region in this plane. 

The correlations between R(D) and R(D*) for the lepton flavor violating operators O e x 
are shown in the right panel of Fig. |4j As these operators do not interfere with the SM one, 
they always increase R(D) and R(D*) with different ratios depending on X. 

Once a future experiment gives precise values of R(D) and R(D*), we may exclude some 
operators depending on the actual experimental values. 
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FIG. 4: Correlations between R(D) and R(D*) in the presence of new physics operators of O t x (left) 
and O'jf' (right) for X = V12, Si 2, T. The black dot in each panel indicates the SM prediction. 
The light blue horizontal and vertical bands are the experimental values. 




R{D) R(D') R(D') 



FIG. 5: Correlations between R(D^) and P T (D^), and R(D*) and P D * in the presence of new 
physics operators O v ^' T (gray horizontal lines), Oy 2 (light gray regions), O v ^ (black dashed lines), 
( ma g en ta curves), OTp (blue regions) and Ojf 1 (white dotted curves). The black dot in each 
panel indicates the SM prediction. The light blue vertical bands show the experimental constraints. 

C. Correlations between decays rates and polarizations 

As shown in Fig. [2j polarizations are also useful observables to identify new physics. Here 
we show correlations between decay rates and polarizations in Fig. [5] The gray horizontal 
lines represent the correlations for Oy^ ,T , the light gray regions (black dashed lines) for 
\ the magenta curves for the scalar operators O e g^ ,T , and the blue regions (white 
dotted curves) for O^S . The operator O v gives the same line as O l v in B — > Dru. The 
light blue vertical bands show the experimental constraints on R(D) and R(D*). 
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We find specific features of the scalar operators in this figure. The polarizations P T (D^*') 
and Pd* are uniquely related to the corresponding decay rates in the presence of scalar 
operators, because the scalar operators O l s and O l S2 contribute only to r + (D(*>) and T(D* L ). 
The definitions of polarizations may be rewritten as 

[r+pW) + r-(D w )] [i - p T (i>W)] = 2r-(,D w ) , (46) 

[T(D* T ) + T(D* L )] (1 -P D *) = T(D* T ) . (47) 

The right-hand sides of these equations are given solely by the SM contributions, and thus 
the polarizations are definitely determined by the corresponding decay rates as seen in Fig. |5j 
These specific correlations are prominent predictions of the scalar type operators, although 
we cannot discriminate O s from O l S2 by using these correlations. 

As for the other operators, it is apparent that the SM operator Oy and Oy do not 
affect the polarizations, while the quark right-handed current O l V2 has no effect on the tau 
polarization in B — > Dtv because the axial vector part does not contribute to this process. 
The operator Oy 2 in B — > D*tv and the tensor operator in both the processes predict 
definite relations between the polarizations and the rates. The operator Oy 2 in B — > D*tv 
and the tensor operator in both the processes have no such specific relations, but the 
covered regions are rather restricted. 

The above correlations that include the polarizations definitely increase the ability to 
restrict possible new physics. We might uniquely identify the new physics operator by these 
correlations in some cases. However their usefulness depends on experimental situations as 
we will see in the next section. 

IV. MODEL-INDEPENDENT ANALYSIS OF NEW PHYSICS 

In this section, we illustrate several possibilities to restrict or identify new physics using 
the observables discussed in the previous section and decay distributions. We suppose that 
R{D) and R(D*) will be measured more precisely in a future super B factory experiment. 
Then we will determine the Wilson coefficient C l x associated with O l x that is assumed to 
be dominant except for the SM contribution and predict other observables. 

Here we consider the following three cases of (R(D), R(D*)): 
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R(D) R{D) R(D) 

FIG. 6: Three examples for the discrimination of new physics contribution: (a) R(D) = 
0.21, R{D*) = 0.29, (b) R{D) = 0.37, R(D*) = 0.28 and (c) R{D) = 0.51, R{D*) = 0.25. 





(a) 


(b) 


(c) 


(R(D),R{D*)) 


(0.21,0.29) 


(0.37,0.28) 


(0.51,0.25) 


X 


s 2 


s 2 


v 2 


T 


Si 


s 2 


v 2 


C T x 


-1.20±i0.18 


-0.81 ±i 0.87 


0.03 ±i 0.40 


0.16±i0.14 


-0.50 ±i 1.08 


0.21 ± i0.56 


0.18±i0.53 


P T (D) 


0.02 


0.44 


0.33 


0.22 


0.60 


0.60 


0.33 




-0.30 


-0.35 


-0.50 


-0.26 


-0.51 


-0.51 


-0.50 


P D * 


0.53 


0.51 


0.45 


0.32 


0.45 


0.45 


0.44 



TABLE I: Predictions for the polarizations in each of three cases. 

(a) (0.21, 0.29) as shown in Fig. |6^a), in which O t S2 is unambiguously identified as the new 
physics operator. 

(b) (0.37,0.28) as shown in Fig. |6^b), in which Og 2 ,Oy 2 , and 0? are the candidates for 
the new physics operator. 

(c) (0.51, 0.25) as shown in Fig. [6|c), in which O r Sl , 0$ 2 , and Oy 2 are possible. 

One of R(D) and R(D*) is chosen to be within the 2a range but the other is allowed to 
deviate more. We note that O^f do not reproduce the assumed sets of R{D) and R(D*). 
Table [T] summarizes the Wilson coefficient and the predicted polarizations for each case. 

In the ), the dominant new physics operator is uniquely determined to be O t S2 . 

The polarizations can be used to confirm that the deviation from the SM comes from the 
operator O t S2 . Furthermore the assumption of one-operator dominance will be tested. 

In the case (b), the dominant operator is 0$ 2 , Oy 2 , or CJ. The predicted values of polar- 
izations vary from operator to operator. We will determine the dominant new physics oper- 
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FIG. 7: Distributions of w in B — > Dtv and B — > D*tu for the case (c), where w = {m 2 B — m 2 M — 
q 2 )/(2mBmM)- The red, magenta, and light gray regions represent the distributions in the cases 



of Og i ,Og 2 and Oy 2 respectively, including the theoretical uncertainties as in Sec. Ill A The w 
distributions in B — > D*tv are almost identical for the three operators. 

ator by measuring P T (D*) for example. Then the one-operator dominance will be checked 
by looking up other polarizations. 

In the case (c), the dominant operator is Og i ,Og 2 , or Oy 2 . The occurrence of Oy 2 is 
distinguished from that of O t s or O t S2 by polarizations. Two scalar operators O r s and 



0$ 2 , however, predict the same values of the polarizations as explained in Sec. Ill C Under 
such a situation, q 2 distributions may discriminate between O t Si and O t S2 . In Fig. [7j we 
present q 2 distributions for T S , Og 2 , and Oy 2 in the case (c). We note that the abscissa is 
w = {m 2 B — m 2 M — q 2 ) / '{2m Bin m) instead of q 2 . The q 2 distribution in B — )■ Dtv turns out 
to be useful for the discrimination in this case. 



V. MODEL ANALYSIS 

In this section, we discuss some new physics models which affect B — > D^rv based on 
the results in Sec. IIHI 

A. 2HDMs and MSSM 

As known well, the charged Higgs boson in two Higgs doublet models (2HDMs) con- 
tributes to B — > D(*'t9 and its effect is enhanced in some cases. In order to forbid flavor 
changing neutral currents (FCNC) at the tree level, a Z 2 symmetry is often imposed in this 
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Type I Type II TypeX TypeY 



£ d cot 2 tan 2 (3 

^ - COt 2 13 1 



cot 2 ft 



TABLE II: Parameters £d u in each type of 2HDMs. 



class of models and it results in four distinct 2HDMs [4TH44] . Their Yukawa terms are 
described as 

Cy = -QlY u H 2 u r - Q L Y d H 2 d R - L L Y e H 2 £ R + h.c. (type I) , (48) 

Cy = -QlY u H 2 u r - Q L Y d Hid R - L L Y t Hi£ R + h.c. (type II) , (49) 

C y = ~QlY u H 2 u r - Q L Y d H 2 d R - L L Y,H^ R + h.c. (type X) , (50) 

Cy = -QlY u H 2 u r - QLY d Hid R - L L Y £ H 2 £ R + h.c. (type Y) , (51) 

where H\ }2 are Higgs doublets defined as 

//; ( K I • //, 'W/;. (52) 




and Vi denotes the vacuum expectation value (VEV) of Hi. The ratio of two VEVs is defined 
as tan/3 = v 2 /v\ and v = \J v \ + v \ = 246GeV. In type I, all masses of quarks and leptons 
are given by the VEV of H 2 . In type II, the up-type quarks obtain their masses from H 2 , 
while the down- type quarks and leptons from Hi. In type X, H 2 and Hi are responsible to 
the quark and lepton masses respectively. The masses of the down-type quarks are given 
by Hi and other fermions acquire their masses from H 2 in type Y. Under this definition v 2 
generates up-quark masses in any type of Yukawa interaction. 

These 2HDMs contain a pair of physical charged Higgs bosons, which contributes to 
B — > D^tu. The relevant Wilson coefficients introduced in Eq. m« are given by 



_ m b m T c _ m c m T c 

in H ± in H ± 
where is the mass of the charged Higgs boson. The parameters £^ and £ u are presented 
in Table [TTj To have a sizable charged Higgs effect, |^,«| should be much larger than unity 
taking the experimental bound on the charged Higgs mass into account. Then the case of 
£ n = 1 or £ rf = — 1 is not acceptable. The case of £ u = — cot 2 P or £ rf = cot 2 P with cot 2 P ^> 1 
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Z 'u 


z d 


Ze 


Type I 


^f-cot/3 


- e u sin/3(l + cot 2 /?) 


- ^ cot /3 + e d sin 0(1 + cot 2 /?) 


^cot/3 


Type II 


^f«cot/3 


e u cos /3(tan /3 + cot /3) 


^f d tan^ e d sin/3(tan/3 + cot/3) 


^tan/3 


TypeX 




-e u sin/3(l + cot 2 /3) 


- ^ cot /3 + e d sin 0(1 + cot 2 (3) 


^ tan 


TypeY 




e u cos /3(tan /3 + cot /3) 


^ Md tanP e d sin/3(tan/3 + cot/3) 


^cot/3 



TABLE III: The matrices Z u ^£ in each type of the 2HDM in the presence of the Z2 breaking terms. 

is unnatural since the top Yukawa interaction becomes nonperturbative. The requirement 
for the top Yukawa interaction to be perturbative results in tan > 0.4 jH]. Therefore, 
£d = tan 2 in 2HDM of type II is suitable and only C T Sl is enhanced. As we have shown in 



Sec. Ill, however, it is difficult to explain the current experimental results by O t s alone. We 
find that this model is disfavored with about 4<r using the combined experimental values in 
Eq. (|2]). This result is consistent with the recent study in Ref. [5]. 

A possible solution within 2HDMs is to violate the Z2 symmetry at the cost of FCNC 
We introduce the following Z 2 breaking terms in the above four models: 

AC Y = -Ql^HxUr - QLe'dHxdR + h.c. (for types I and X) , (54) 
A£ Y = —Ql£' u Hiu r — Ql^J^^R + n - c - (f° r tyP es II an d Y) , (55) 

where e' u d are 3x3 matrices that control FCNC and the quark fields are those in the weak 
basis. Writing the Yukawa terms Cy + ACy in terms of mass eigenstates, we obtain the 
following physical charged Higgs interaction terms: 

C H ± = [u R Z\y CKm dL + u L V C KuZdd R + v L Z t l R ) H + + h.c. , (56) 



where Vckm is Cabibbo-Kobayashi-Maskawa (CKM) matrix. Table III shows the expressions 
of Z u ^,e, where M u ^t denote the diagonal up-type quark, down-type quark, and lepton mass 
matrices, and e U) d represent matrices e' u d in the quark mass basis. 

The FCNC in the down-quark sector is strongly constrained, so that is negligible in 
the present analysis. On the other hand, constraints on the FCNC in the up quark sector 
are rather weak. Recently the 2HDM of type II that allows FCNC in the up quark sector is 



studied to explain B — > D^ru jlQ]. Table III implies that the operator O t S2 in types II and 



X might be significant for large tan /3. Then the corresponding Wilson coefficient is given 
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by 

Cs^^-^ieir sm/3 tan/3. (57) 



As seen in Fig. [3] the current experimental results are described by the 2HDM of the type II 
or X with FCNC provided that |e* c | ~ 1. If this is the case, we expect sizable deviations in 
polarizations P T (D^) and Pp* from the SM as designated by the magenta curves in Fig. [5} 
The charged Higgs effects on B — > D^tv in the minimal supersymmetric standard model 
(MSSM) are the same as those in the 2HDM of type II at the tree level. Loop corrections 



induce non-holomorphic terms e U) d in Eq. (55) [19, H5]. However it seems difficult to enhance 
e„ c to be 0(1). Thus the sufficient enhancement of O t s is unlikely in MSSM. 



B. MSSM with R-parity violation 

The ^-parity violating (RPV) MSSM [H] may also affect B D^rv HMH]. We 
consider the following superpotential: 

Wrpv = ^ijkLiLjEl + \' i , jk L i QjDl , (58) 

where A^ and \'^ k are RPV couplings and i,j, k are generation indices. Apart from the 
charged Higgs contribution, there are two kinds of diagrams which contribute to B — > D^ru, 
that is, the slepton and down squark exchanging diagrams. The corresponding effective 
Lagrangian is written as 



4T = -5> 

j,k=l 



3 

T 2fe 



A ij3*jk3 - , - i . A i3j A 3kj - , c \ /-c\i , 

Y—c L b R T R v % L + — ~^c L r C ) R {v c ) R b L 

mi mi 

l L d R 



(59) 



where mp (m$ ) is the mass of the slepton (down squark) for the j-th generation and Vy 
is the component of CKM matrix. Here we assume that the slepton and down squark mass 
matrices are diagonal in the super-CKM basis for simplicity. Using Fierz identity the second 



term in Eq. ( 59 ) is rewritten as 



c L (r% (9% b L = \c Ll »b L f LlfM ul (60) 
Then, the corresponding coefficients, defined in Eq. Q, are expressed as 

18 



Ill 



On the 



We find that \C Vi \ > 0.08 and 



where l\ — e, I2 — /i, and h = t represent the flavor of the neutrino. 

The case that Cg is dominant contribution is disfavored as discussed in Sec. 
other hand C v has an allowed region as shown in Fig. 
|CyH > 0.42 (90% CL.) are required to explain the current experimental results. However 
the same products of RPV couplings as in Cy 's also contribute to B — > X s vv (50]. The 
upper limit of the branching ratio, B(B — > X s uu) < 6.4 x 10~ 4 (90% CL.) [51], puts a strong 
constraint on RPV couplings. It turns out that (CyJ < 0.03 is imposed, which contradicts 
the above requirement for B — > D^'tp. Thus MSSM with RPV is not likely to be consistent 
with both B — > D^'tv and B — >■ X s vv at the same time. 



C. Leptoquark models 

There are ten possible leptoquark models that respect the symmetry of the SM |52j. 
Among them, the following leptoquark model is interesting in the sense that the tensor 
operator is generated [53] . 

Clq = (XijQ^R + , (62) 

where Slq is the scalar leptoquark with SU(3) x SU(2) x Z7(l)y quantum number being 
(3, 2, 7/6). The effective Lagrangian is represented as 

- 4? = ^ (f R b L ) (c R vl) = -^k (c R b L f R vl + \c R a^b L r R a, v v\\ . (63) 

Both 0% and appear at the same time, however the Wilson coefficients are related to 
each other: 

* _ A nk _ - 1 A 3 3 A' 2 , ^ n R fXh^\ ^500GeV\ 2 



C'* = 4C£ = — - -0-6 . (64) 

32 2V2G F V cb 2m| LQ V 0.4 A ™s LQ J 

Thus, a similar analysis as the case of one dominant operator can be made. 

The bottom right panel in Fig. [3] shows the constraint on the above Wilson coefficient 
Cg 2 (= 4CJ). We find that the present experimental results are explained with Cg 2 ~ ±i 0.6 
or [Cg^l ~ 0.6. The relevant mass scale of the leptoquark could be ~ 500 GeV as seen in 



Eq. (64). In Fig. |8j we show the correlations between decay rates and polarizations in the 
above leptoquark model. The green regions represent the correlations, R(D^)-P T (D^) and 
R(D*)-Pd*, as in Fig. [5} Taking into account both the experimental constraints on R(D) 
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FIG. 8: Correlations between R(D^) and P T (D^), and R(D*) and Pd* in the leptoquark model 



of Eq. (62) are represented by green regions. Yellow regions indicate the constraints from both the 



present experimental bounds on R(D) and R{D*). The black dot in each panel stands for the SM 
prediction. 



and R(D*) at the same time, we find the present allowed regions at 99% C.L. as shown by 
the yellow regions. 



VI. CONCLUSIONS 

We have studied the exclusive semi-tauonic B decays, B — > D^*'ri>, in the model- 
independent manner based on the effective Lagrangian including all the possible four-Fermi 
operators. It has turned out that the current experimental values of R(D) and R(D*) are 
not explained by the operator O l Sl nor O e £ alone, while the other operators Og 2 , O l Vi , O l y 2 , 
and O l T reasonably work under the assumption of one-operator dominance. More precise 
data that will be given in a future super B factory experiment will allow us to identify the 
relevant new physics operator among these operators if the deviation from the SM persists. 
We have pointed out that correlations among observables including the longitudinal tau 
polarizations, the D* polarization and q 2 distributions are useful in distinguishing among 
possible new physics. This observation for the most general effective Lagrangian in Eq. Q 
confirms and extends results in the literature for limited new physics cases [25j 123 I2H1 GFJ] . 

Furthermore, we have studied several interesting models that may describe the deviations 
of R{D) and R(D*) from the SM based on our model-independent analysis. In 2HDMs 
without tree-level FCNC and MSSM, only T S could be enhanced and thus are unlikely to 
explain the deviations, while T S2 generated in 2HDMs of type II and type X with FCNC 
might be sizable depending on the magnitude of the t — > c FCNC parameter. The parameter 
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region of MSSM with i?-parity violation that explains the deviations in R(D) and R(D*) 
has turned out to be inconsistent with the experimental bound on B — > X s uu. The scalar 
leptoquark model that simultaneously induces O l s and O l T describes the experimental data. 

In conclusion, B — > D^tv is a powerful tool to explore new physics in the charged 
current. It is important to study combinations of observables including particle polarizations 
and q 2 distributions as well as the decay rates in order to clarify possible new physics 
in a model-independent manner. The model-independent analysis gives a firm footing in 
examining models of new physics. 
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Appendix A: Leptonic amplitudes 

In this appendix, we summarize the expressions of the leptonic amplitudes. The vector 
type leptonic amplitudes are given by 

T VSmd T , (Al) 

Lq = 2m T vcos9 T , (A2) 

L+ = -2m T v , (A3) 

L± = V2^v(l±cos9 T ) , (A4) 

L = -2^vsm9 T , (A5) 

L~ = , (A6) 

where v = \J\ — m^/q 2 . The scalar type leptonic amplitudes are written as 

L + = -2^v , (A7) 

L~ = . (A8) 



21 



The tensor type leptonic amplitudes are 

^aa = , (A9) 

L o± = ~ L to = -V2y^vsm9 T , (A10) 

L%_ = -L± + = L+ = -L+ = 2y^vcos9 T , (All) 

L+ s = -L + s± = ^V2^/fvs\n6 T , (A12) 

Lq± = -L ±0 = ^V2m T v(l ±cos8 T ) , (A13) 

L^__ = —LZ + = L~ Q = —Lq s = —2m T v sin 9 T , (A14) 

L~ s = -LJ ± = -V2m T v(l±cos9 T ) . (A15) 

The subscript I that represents the neutrino flavor is omitted for brevity. 
Appendix B: Form factors and hadronic amplitudes 

Here, we show the definitions of the relevant form factors and the hadronic amplitudes 
in B — > D^tu. We also discuss how to evaluate the scalar and tensor form factors by using 
the equations of motion. 

1. Vector and axial vector operators 

We define the form factors of the vector and axial vector operators in B — > D^*'ru as, 

(D(p D )\crfb\B(p B )) = ^^[h + (w)(v + v'r + h4w)(v-v'r} , (Bl) 

(D*(p D .,e)\crfb\B(p B )) = iy/m B m D * h v (w)e^ pa etv' p v a , (B2) 
(D*(p D .,e)\crfy s b\B(p B )) = y/m B m D * [h Al (w)(w + l)e*" 

-(e* • v) (h Aa (wX + h A3 (w)v'») ] , (B3) 

where v = Pb/^b.v' = Pm/ttim, and w = v ■ v' in B — > Mrv. The polarization vector of 
D* is indicated as e M . 



Using these definitions, the hadronic amplitudes defined in Eqs. (18) and (19) are repre 
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sented as 



H v u ± = H k,± = 



Vi,0 



^Jvfi - 1 

#v 2 ,o = m B v^^^ [(1 + r)h+(w) - (1 - r)/i_(w)] 



(B4) 
(B5) 



H S V2 . S = m B Vr^== [(1 -r)(w + l)h + (w) -(l + r)(w- 1)M«0] , (B6) 
q z {w) 



in B — >■ £)ri/ and those in 5 — > D*tu are given by 



m B y/r (w + l)h Al (w) =F Vw 2 - l/ly(w) 



mB\fr 



w + 1 



(w - r)h Al (w) + (w- l)(rh A2 (w) + /u 3 (w))] , 



(B7) 



(B8) 



~ H V 2 ,s 



= m B ^/r—== [-(w + l)h Al (w) + (1 - rw)h M (w) + (w - r)h Aa (w)} , (B9) 
others = 0. (BIO) 



In the heavy quark limit, we obtain h + (w) = h Al3 (w) = hv{w) = and h_(w) = 

h A2 (w) = where £(w) is Isgur-Wise function [38]. In order to evaluate these form factors 
by using dispersion relations and the heavy quark effective theory, it is convenient to redefine 
the form factors as follows [36J: 

1 



V x (w) 
Sx(w) 

Ri(w) 
R 3 (w) 



1 + r v ' 
1 + r w — 1 

1-rw+l 



h+(w) - 

h+(w) - 

h Al (w) , 
h v (w) 
h Al (w) ' 

h A . A {w) +rh M {w) 

h Al (w) 
h A3 (w) - rh M (w) 

h Al (w) 



h-(w) 



(Bll) 

(B12) 
(B13) 
(B14) 

(B15) 

(B16) 
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2. Scalar and pseudo-scalar operators 



For the scalar and pseudo-scalar operators, the form factors are defined as 



(D(p D )\cb\B(p B )) = ^m B m D (w + l)h s (w) , 
(D*(p D .,e)\cry 5 b\B(p B )) = -^/m B m D » (e* ■ v)h P (w) 



and then the hadronic amplitudes in Eqs. (20) and (21) are represented as 



H k = HL = m Bv ^{w + l)h s (w) , 



in B — > Dtu, and 



H Si = H s 2 =0 



H v Si = -H U S2 = -m B V^Vw 2 -lh P (w 

in B -> D*tv. 

Using the quark equations of motion we obtain 

2<9 M (c7 M 6) = (m b - m c ) cb , 
{dp (c7 M 7 5 6) = - (m b + m c ) C7 5 fe, 

which lead to the following relations among form factors: 



h s (w) 



h+(w) 
SAw)- 



r w 



O 



A 2 



h-(w) + 



m. 



hp(w) 



— |— [{w + l)h Al {w) + (rw - l)h A2 {w) - (w - r)h A . A {w)} + O ( — j 
1+r ~ \ m Qj 



A 1 (w) 
1 + r 



w + 1 



2r 



2r 



°(-) 

\ m Qj 



(B17) 
(B18) 



(B19) 



(B20) 
(B21) 



(B22) 
(B23) 



(B24) 



(B25) 



where we have used m B {M) — m b{c) + A + O (A 2 /m 6 ( c )). The absence of l/mq corrections 



in Eq. (B24) is confirmed by the heavy quark expansion without resorting to the equations 



of motion. Thus we employ hs(w) = Si(w) = [1 + A(w)] Vi(w), where A(w) involves the 



1/mg corrections. On the other hand, there exist unknown I/tuq corrections in Eq. (B25). 



We ignore them, but use the axial vector form factors Ai(w) and R2${w) including I/thq 
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corrections as described in Sec. II C The I/ttlq corrections in R,2(w) are necessary to evaluate 
the rate of B — > D*£u, because the experimental determination of the slope parameter p\ 
assumes these corrections EZ|. As for R 3 (w), since it does not contribute to B — > D*£u and 



no experimental data are available, we use the theoretical estimation by the heavy quark 
effective theory with ±100% error in I/tuq corrections. 



3. Tensor operator 



The form factors of the tensor operator are defined as 



(D(p D )\ccx^b\B{p B )) = -i^m B m D h T {w){v fl v' u - v v v'») , (B26) 
(D*(p D *,e)\ca^b\B(p B )) = -^m B m D *e^ \h Tl (w)e* p (v + v% + h T2 (w)e* p (v - v') a 



+h T3 (w)(e* -v)(v + v') p (v - v') a 



(B27) 



where a^ v = {i/2) [7^, 7"] and e 0123 = —1. The amplitudes corresponding to the operator 
ccr MV 756 are given by the following identity, 



^"75 



(B28) 



With the above form factors, the hadronic amplitudes in Eq. (22) are represented as 



HQ S (q 2 ) = —m B y/r\/w 2 — 1 hr(w) 
±Hl(q 2 ) 



(B29) 



±m B \fr 



1 — rw ± r\Jw 2 — 1 



x 

lJ2 



h Tl (w) + h T2 {w) + (w ± \/w 2 - l){h Tl {w) - hr 2 {w)) 



(B30) 



HUq 1 ) = H° s (q 



-m B 



r [{w + l)h Tl {w) + (w- l)h T2 {w) + 2(w 2 - l)h Ta {w)] , (B31) 



AA' 

others 



a' a vy ; ) 



(B32) 
(B33) 



As in the case of the scalar type operators, the equations of motion lead to 



d.(ca^b) = -{m b + m c )cYb-{id v c)b + c{id v b) 



(B34) 
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and we obtain 



h T (w) = h+(w) - \^-h-(w) + (—) , (B35) 



1 — r \ m Q J 

h Tl (w) = -^[{l-r) 2 {w + l)h Al {w)-{l + r) 2 {w-l)h v {w)]+o(—),{B3^ 
2q 2 (w) \ m Q/ 

h T2 (w) = (1 ~ j! 2 )(w + 1} [h Al (w) - h v (w)] + (—) , (B37) 
2q 2 (w) V m Q/ 

hr s {w) = ~ 2 (i + l r )qi(w) [ 2r( ^ w + l ) hA ^ w ) ~ ( r/ ^ 2 ( w ) - ^M) 

-(l+f) 2 ^Hl +0^— V (B38) 



We find hriw) = h^iw) = £(w) and h,T 2 (w) = hr 3 {w) = in the heavy quark limit. 
Similarly to the case of the scalar type operators, we ignore the unknown l/mq corrections 
in the right-hand sides, and employ the vector and axial vector form factors with the 1/tuq 
corrections. 

Appendix C: Decay rates 

The differential decay rate B — > Mrv is represented as 

d^t = ^"E\My»tf,casO T )\ 2 d* 3 , (CI) 

where the three-body phase space d$3 is given by 

rf$ 3 = ^\ Q - (l-^dq'dcosOr, (C2) 



2567T 3 m 2 B \ q 2 / 

and Q± = {tub ± itlm) 2 — <f- Several decay rates used in the main text are obtained by 



integrating Eq. (CI) over q 2 and cos# T . For notational simplicity, we define the following 
quantities: 

T ± {D) = T ±i/2 } T ± m= rtf, (C3) 



A M =±1,0 
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